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Abstract 

We study the robustness of the determination of the neutrino masses and mixing from the 
analysis of atmospheric and K2K data under the presence of different forms of phenomenologically 
allowed new physics in the Vp-v T sector. We focus on vector and tensor- like new physics interactions 
which allow us to treat, in a model independent way, effects due to the violation of the equivalence 
principle, violations of the Lorentz invariance both CPT conserving and CPT violating, non- 
universal couplings to a torsion field and non-standard neutrino interactions with matter. We 
perform a global analysis of the full atmospheric data from SKI together with long baseline K2K 
data in the presence of v„ — > v T transitions driven by neutrino masses and mixing together with 
sub-dominant effects due to these forms of new physics. We show that within the present degree 
of experimental precision, the extracted values of masses and mixing are robust under those effects 
and we derive the upper bounds on the possible strength of these new interactions in the v^—v T 
sector. 
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I. INTRODUCTION 



Neutrino oscillations are entering a new era in which the observations from underground 
experiments obtained with neutrino beams provided to us by Nature - either from the Sun 
or from the interactions of cosmic rays in the upper atmosphere - are being confirmed by 
experiments using "man-made" neutrinos from accelerators and nuclear reactors Q]. 

For atmospheric neutrinos, Super-Kamiokande (SK) high statistics data 0,0] established 
beyond doubt that the observed deficit in the //-like atmospheric events is due to the neutri- 
nos arriving in the detector at large zenith angles, and it is best explained by oscillations. 
This evidence was also confirmed by other atmospheric experiments such as MACRO 0] 
and Soudan 2 0. 

The KEK to Kamioka long-baseline neutrino oscillation experiment (K2K) uses an 
accelerator-produced neutrino beam mostly consisting of with a mean energy of 1.3 GeV 
and a neutrino flight distance of 250 km to probe the same oscillations that were explored 



with atmospheric neutrinos. Their results [oj show that both the number of observed neu- 
trino events and the observed energy spectrum are consistent with neutrino oscillations, with 
oscillation parameters consistent with the ones suggested by atmospheric neutrinos. 

Oscillations are not the only possible mechanism for atmospheric — > v T flavour transi- 
tions. They can also be generated by a variety of forms of nonstandard neutrino interactions 
or properties. In general these alternative mechanisms share a common feature: they require 
the existence of an interaction (other than the neutrino mass terms) that can mix neutrino 



flavours J7(. Among others this effect can arise due to violations of t 



(VEP) non-standard neutrino interactions with matter neutrino couplings 

to space-time torsion fields |12|, violations of Lorentz invariance (VLI) 0, Q and of CPT 



re equivalence principle 



symmetry |15l Il6l|. From the point of view of neutrino oscillation phenomenology, the most 
relevant feature of these scenarios is that, in general, they imply a departure from the E^ 1 
energy dependence of the oscillation wavelength jl7| . 

Prior to the highest-statistics SK data, some of these scenarios could provide a good 
description - alternative to Am 2 neutrino oscillations - of the atmospheric neutrino phe- 
nomenology 0, 0] . However, with more precise data, and in particular with the expan- 
sion of the energy range covered by atmospheric neutrino data due to the inclusion of the 
upward-going muons, these alternative scenarios became disfavoured as leading mechanism 



to explain the observations 20j, |2jJ, |22| . The results from K2K experiment further singled 
out oscillations as the dominant mechanism of <-> v T transitions j^. 1 

The question arises, however, to what point the possible presence of these forms of new 



1 Recently |3 SK collaboration has presented a reanalysis of the SKI data in terms of the reconstructed 
L/E which allowed them to slightly improve the discrimination between oscillations and alternative mecha- 
nisms. Unfortunately to reproduce such analysis for the subdominant effects discussed here is not possible 
outside the collaboration. 
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physics, even if sub-dominant, may affect the derived ranges of masses and mixing from the 
oscillation analysis of the atmospheric and K2K data. Or in other words, to what level our 
present determination of the neutrino masses and mixing is robust under the presence of 
phenomenologically allowed new physics effects. 

In this paper we address this question by performing a global analysis of the atmospheric 
and K2K data with — > v T transitions driven by neutrino masses and mixing in the pres- 
ence of some generic forms of new physics. In particular we consider new physics interactions 
which are vector-like, or tensor-like (scalar interactions cannot be distinguish from oscilla- 
tions). This allow us to treat, in a model independent way, effects due to the violation of the 
equivalence principle, violations of the Lorentz invariance both CPT conserving and CPT 
violating, non-universal couplings to a torsion field and non-standard neutrino interactions 
in matter. In Sec. |Il] we present the formalism adopted and the data set used. In Sec. IIHI 
we show the results of our analysis. Conclusions are given in Sec. IIVI The technical details 
of our new statistical analysis of the atmospheric data are described in the appendix. 



II. FORMALISM 



In what follows we consider some new physics (NP) scenarios which induce new sources 
of lepton flavour mixing in addition to the "standard" Am 2 oscillations (Am 2 -OSC). We 
concentrate on flavour mixing mechanisms for which the propagation of neutrinos (+) and 
antineutrinos (— ) is governed by the following Hamiltonian jl6| |: 

Am\ T f-1 o\ ± i A5 n E n ( -I ()\ TT , 

where Am 2 is the mass-squared difference between the two neutrino mass eigenstates, 
accounts for a possible relative sign of the NP effects between neutrinos and antineutrinos 
and A5 n parametrizes the size of the NP terms. The matrices U# and "U^ nt ± Vn are given by: 

tt = [ cos # sind \ tt ( cos £n sin £ n e ±ir iA 

6 \-sm6 cosOj ' Cn ' ±r?n ^-sin£ n e^« cos£ n J ' 

where we have also accounted for possible no n- vanishing relative phases r] n . For concreteness 
we will focus on NP effects which are induced by tensor-like and vector-like interactions. 

We denote by tensor-like interactions those with n — 1 leading to a contribution to the 
oscillation wavelength which grows linearly with the neutrino energy. For example, Eq. (JT|) 
can describe the evolution of and u T 's of different masses in the presence of violation of 
the equivalence principle (VEP) due to non universal coupling of the neutrinos, 71 7^ 72 
(i/x and u 2 being related to and v T by a rotation 9q), to the local gravitational potential 
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j. 2 Phenomenology of neutrino oscillations induced or modified by VEP has been 
widely studied in the literature [2^ . 
In this case 

A5 1 = 2|0|( 7l -7 2 ) = 2|0|A 7 , £i = G , °l = °i- (3) 

For constant potential 0, this mechanism is phenomenologically equivalent to the breakdown 
of Lorentz invariance induced by different asymptotic values of the velocity of the neutrinos, 



v\ 7^ v 2 , with Vi and v 2 being related to and v T by a rotation 9 V 
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. In this case 



A5i = (v! - v 2 ) 5v, £i = 9 v , = a x . (4) 

We denote by vector-like interactions those with n = which induce an energy inde- 
pendent contribution to the oscillation wavelength. This may arise, for instance, from a 
non-universal coupling of the neutrinos, k\ 7^ k 2 (u^ and v 2 being related to the and v T 
by a rotation 9q), to a space-time torsion field Q [12|, so 

A5 = Q(h-k 2 ) = QSk, £ = 9 Q , ct+ = ct . (5) 

Violation of CPT due to Lorentz- viola ting effects also lead to an energy independent con- 
tribution to the oscillation wavelength \lA with 

A5 = b 1 -b 2 = 5b, £ = 0<zpt , o"cf = -&o ( 6 ) 

where 6j are the eigenvalues of the Lorentz violating CPT-odd operator ^fb^f^yU^ and 
9 V is the rotation angle between the corresponding neutrino eigenstates and the flavour 
eigenstates fl^ ]. 

In all these scenarios, if the NP strength is constant along the neutrino trajectory, the 
expression of P v -* Vll takes the form 

P v ^ = 1 - P v ^ Vt = 1 - sin 2 29 sin 2 Kj . (7) 

where the correction to the Am 2 -OSC wavelength, TZ, and to the global mixing angle, 0, 
verify 

K cos 26 = cos 29 + ^2 R n cos 2£ n , (8) 

n 

K sin 26 = |sin 29 + ^ R n sin 2£„ e iVn | , (9) 



2 VEP for massive neutrinos due to quantum effects discussed in Ref. 01 can a l so ^> e parametrized as 
Eq. (JTJ with n = 2. 
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with R n being the ratio between Am 2 -induced and the NP-induced contributions to the 
oscillation wavelength 

, A5 n E n AE , . 

*" = ff »^A^- (10) 
For Py^Vf, the same expressions hold with the exchange cr+ — > cr~ and r\ n — > — rj n . 

For the sake of simplicity, in what follows we concentrate in scenarios with one NP source 
characterized by a unique n. In this case 

sin 2 26 = — (sin 2 29 + R 2 n sin 2 2£ n + 2R n sin 26 sin2£ n costfo) , (11) 
7Z 2 



K= y/l + B* + 2i^ (cos 20 cos 2£ n + sin 20 sin 2^008^) . (12) 

In Fig. ^ we illustrate the effect of the presence of the NP in the atmospheric neutrino 
events distributions for Am 2 -0SC plus sub-dominant CPT-even tensor-like and vector-like 
NP effects, for some characteristic values of the NP-parameters. In both cases R n is a 
growing function of E and the NP effects become relevant in the higher energy samples, in 
particular for upward going muons. 

In order to properly define the intervals of variation of the five parameters Am 2 , 0, A5 n , 
6n> i]n, we can take advantage of the symmetries of the Hamiltonian (see also Ref. [2q for a 
very similar problem) . For a given value of <r+ , from the expressions (Q) and (J2J) we see that 
the Hamiltonian is invariant under the following transformations: 



• £n — > £n + 7T, 

• Vn-*Vn + 27T, 

• Am 2 ^ -Am 2 and 9^9 + n/2, 

• A5 n -> -A5„ and £„ ^ £„ + vr/2, 

• £n -> and Vn -> + 7T- 

Furthermore, the relevant survival probabilities P y _>„ and P Pfi _>p M are not affected by a 
change in the overall sign of the Hamiltonian, as well as change in the global phase of its 
non-diagonal components. Therefore, we also have: 

• 9 -> 9 + vr/2 and £ n -> £ n + vr/2, 

• 9 —> -9 and -> 

• T)n ~> -?7n- 
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FIG. 1: Zenith-angle distributions (normalized to the no-oscillation prediction) for the Super- 
Kamiokande //-like events. The full line gives the distribution for the best fit of Am 2 -OSC, Am 2 = 
2.2 x 10~ 3 eV 2 and sin 2 9 = 0.5. The dashed and dotted lines give the distributions for Am 2 - 
OSC+NP scenarios for n = 1 and n = with A$i = 2.0 x 10" 24 and A5 = 4.2 x 10" 23 GeV 
respectively. In both cases 77 = £ = and the oscillation parameters have been set to their best fit 
values. 



The above set of symmetries allows us to define the ranges of variation of the five parameters 
as follows: 

(a) Am 2 >0, (c)0<#<7r/2, 
(&)A5 n >0, (d) 0<£n<7r/4, (13) 

(e) < i] n < 7r . 

Thus in the general case we cover the mixing parameter space by using, for instance, < 
sin 2 6 < 1 and < sin 2 2£ n < I. 

For the case of real relative phase, rj n G {0, 7r}, one can absorb the two values of r) n into 
the sign of £ n . In this case we drop (e) and replace (d) by: 

(d') -7r/4<£n<7r/4 (14) 

and use instead — 1 < sin2£ n < 1. 
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Finally we notice that the above derivation is valid for a given sign of cr+. Keeping the 
convention of Am 2 > and A5 n > the survival probability for the opposite sign can be 
obtained by the exchange 

sin 2 9 — > 1 — sin 2 9 and i] n ~^ir — r] n . (15) 



In addition, we also consider the special case of vector-like NP due to non-standard 
nentrinc-matter interactions (NSI) QQ. In this case the effective Lagrangian describing 
the evolution of the v ii —v T system can be written as 



iai2 

where Nf(r) is the number density of the fermion / along the path r of the neutrinos 
propagating in the Earth, and e± and e'± parametrize the deviation from standard neutrino 
interactions: \f2GpNf(r) e + is the forward scattering amplitude of the FC process u^ + f — > 
v T + f, and V2G F N f (r) e' + is the difference between the v T + f and the v^ + f elastic forward 
scattering amplitudes. The corresponding amplitudes for antineutrinos are given by e_ and 
e'-. For simplicity we assume that the interactions for neutrinos and antineutrinos are the 
same, which implies e' + = e'_ = e' and e + = e*_ = s. Thus the NSI Hamiltonian contains 
three real parameters, which can be chosen to be s', \e\ and arg(e). NSI and their interplay 
with the oscillations have also been studied in different contexts: among others, in relation to 
snpexnova nentrinos 0, to the solar nentrino problem & to the LSND resnlts oscillation 
results 
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, to neutrinoless double beta decay |3Cj, and to present and future laboratory 

neutrinos 

Formally Eq. (JTfij) can be seen as a special case of Eq. (|TJ) with n = 0, a Q = —Cq, and 
A5 = 2V2 G F N f (r) T = 4.58 x KT 22 (2 - Y p ) P jp^ T GeV , 



3g/cm 3 



cos(20 = ^, sin(20 = ^, r ? = arg(e), (17) 



e' 2 



with JF = y \e\ 2 + — . 

Technically the main difference is that NSI only affect the evolution of neutrinos in the 
Earth, and their strength changes along the neutrino trajectory. Consequently the flavour 
transition probability cannot be simply read from Eq. (J7J) and its evaluation requires the 
numerical solution of the neutrino evolution in the Earth matter. In our calculations we use 
PREM (33! for the Earth's density profile and a chemical composition with proton/nucleon 
ratio Y p = 0.497 in the mantle and 0.468 in the core. In what follows for the sake of 
concreteness we set our normalization on these parameters by considering that the relevant 
neutrino interaction in the Earth occurs only with down-type quarks. 
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Concerning the data samples used in the analysis, for atmospheric neutrinos we include 
in our analysis all the contained events as well as the the upward-going neutrino-induced 
muon fluxes from the latest 1489 SK data set Q . This amounts for a total of 55 data points. 
Details of our new statistical analysis, introduced here for the first time, are presented in 
the appendix. 

For the analysis of K2K we include the data on the normalization and shape of the 
spectrum of single-ring /i-like events as a function of the reconstructed neutrino energy 0|. 
We bin the data in five 0.5 GeV bins with < E rec < 2.5 GeV plus one bin containing all 
events above 2.5 GeV. Details of our analysis of the K2K data were presented in Ref. j^l and 
will not be repeated here. Let us just comment that together with statistical uncertainties 
we also account for the systematic uncertainties associated with the determination of the 
neutrino energy spectrum in the near detector, the model dependence of the amount of nQE 
contamination, the near/far extrapolation and the overall flux normalization. 



III. RESULTS AND DISCUSSION 



We now describe the results of our x 2 analysis of the standard Am 2 -OSC+NP scenarios. 
As discussed in the previous section the analysis contains five parameters: Am 2 , 9, A5 n , £ n 
and 7] n (or e', \e\ and arg(e) for NSI). Our results are summarized in Figs. El El IH and 
where we show different projections of the allowed 5-dimensional parameter space. 

In Figs. El and El we plot two-dimensional projections of the allowed parameter region for 
the analysis of atmospheric and K2K data in presence of — > v T oscillations and different 
NP effects parametrized in the form Eq. (0). The corresponding results for the case of NSI 
are presented in Fig. El The regions in each panel are obtained after marginalization of x 2 
with respect to the three undisplayed parameters and they are defined at 90%, 95%, 99% 
and 3(7 CL for 2 d.o.f. (A X 2 = 4.61, 5.99, 9.21 and 11.83, respectively). 

The left panels in Figs. |2]andElshow the projection of the allowed region on the oscillation 
parameters Am 2 -sin 2 6 plane. The best fit point is marked with a star. For the sake of 
comparison we also show the lines corresponding to the contours of the allowed regions for 
pure Am 2 -OSC and mark with a triangle the position of the best fit point. The results are 
shown for the chosen relative sign cr+ = +1. For er+ = — 1 the corresponding region would 
be obtained by sin 2 6 — > 1 — sin 2 6. 

The regions on the right panels of Fig. El and El show the allowed values for the parameters 
characterizing the strength and mixing of the NP. The full regions correspond to arbitrary 
values of the relative phase rj n (or equivalently to complex e parameter for the NSI case) 
while the lines show the results for real relative phase rj n e {0, 7r} (which for NSI corresponds 
to £ real and either positive or negative, respectively). For this second case we show the 
allowed region for —1 < sin2£ n < 1 where for cx+ = +1 the sector with —1 < sin2£ n < 
and < sin2£„, < 1 correspond to i] n = n and r\ n = 0, respectively, while the opposite 
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FIG. 2: Allowed parameter regions for the analysis of atmospheric and K2K data in presence of 
v — ► v T oscillations and different NP effects as labeled in the figure. Each panel shows a two- 
dimensional projection of the allowed five-dimensional region after marginalization with respect 
to the three undisplayed parameters. The different contours correspond to the two-dimensional 
allowed regions at 90%, 95%, 99% and 3<r CL. The filled areas in the left panels show the projected 
two-dimensional allowed region on the oscillation parameters Am 2 -sin 2 6 plane. The best fit point 
is marked with a star. For the sake of comparison we also show the lines corresponding to the 
contours in the absence of new physics and mark with a triangle the position of the best fit point. 
The results are shown for the chosen relative sign <r+ = +1; for <r+ = — 1 the corresponding region 
would be obtained by sin 2 — > 1 — sin 2 9. The regions on the right panels show the allowed values 
for the parameters characterizing the strength and mixing of the NP. The full regions corresponds 
to arbitrary values of the phase rj n while the lines correspond to the case r) n E {0, it}. 
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FIG. 3: Same as Fig. for the case of Am 2 -OSC+NSI. 



holds for <r+ = — 1. As discussed in the previous section, for the case of arbitrary phase r] n 
the full mixing parameter space can be covered by < £ n < 7r/4, which translates into the 
symmetry of the allowed region around £ n = 0. 

Several comments are in order. First, from the figures we see that the best fit point for 
the global Am 2 -OSC+NP scenarios is always very near the best fit point of pure Am 2 -OSC 

Am^ est = 2.2 x 1(T 3 eV 2 sin 2 # best = 0.5 . (18) 

In other words, the data does not show any evidence of presence of NP even as a sub- 
dominant effect. Second, in agreement with SK analysis [2], we find that with the inclusion 
of the three-dimensional atmospheric fluxes and improved cross sections as well as with 
the reanalyzed data points from SK, the best fit point and corresponding allowed regions 
from the atmospheric+K2K neutrino analysis is shifted to slightly lower values of Am 2 
compared to our previous analysis corresponding to the same data set j^. Third, the 
figures illustrate the robustness of the allowed ranges of mass and mixing derived from the 
analysis of atmospheric and K2K data under the presence of these generic NP effects. Fourth 
the analysis allow us to derive well-defined upper bounds on the NP strength. 

From Fig. Elwe see that the bounds on the NP strength parameter A5 n tightens for larger 
values of £ n , being this effect stronger for NP effects leading to sub-dominant oscillations 
with stronger energy dependence. In particular, for n = 1 the bound on A5 n for sin 2 2£ n = 1 
is a factor ~ 50 stronger than that for £ n = 0, while for n = the variation of the bound 
on A5 n with £ n is at most a factor ~ 3. This behaviour can be qualitatively understood by 
studying the modification of the oscillation probability at the best fit point of oscillations, 
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best 



1, due to NP effects: 



AP 





1 + 2R n>b sin2£ n cosr? n + R 2 n b sin 2 2£, 



)Am 2 



1 + 2P n>b sin 2£ n cosr? n + 



■) 



(19) 



sin 2 U b Jl + 2R nfi sin 2£ n cos r/„ + R 2 n b 



1 — sin 2 b 



where (fib = 2.8 (L/10 3 km) (GeV/E) is the Am 2 oscillation phase at the best fit point and 
R nyb = 0.91 x 10 21 (A5„/GeV^ n ) (P/GeV) n+1 is the ratio of NP to the standard oscillation 
contributions evaluated at the best fit point of oscillations. 

From Eq. (fT9"|) we find that as long as (fib is small the dependence of AP on R n ^ (and 
consequently on A6 n ) is stronger for larger values of | sin2£J, which explains the tightening 
of the bound on AS n . This behaviour was found in Ref. [20| for the case with n = 1. 

However, it is worth noticing that the characteristic value of (fib for which NP effects are 
relevant depends on n since as n increases the effect becomes important only for higher 
values of E (see Fig. HJ. As a consequence, the characteristic (fib for n = is larger than for 
n = 1. Numerical inspection of Eq. ()19j) also shows that the variation of the dependence of 
AP on R Ut b with sin 2£ n decreases as (fib increases. This explains the milder dependence of 
the bound on A5 n with the mixing angle sin 2£ n for n = as compared with the n = 1 case. 
The figure also illustrates that imposing that the Hamiltonian is real does not substantially 
affect these conclusions. 

The same arguments apply to the results for NSI in Fig. El In particular one sees that, 
as expected, the results for NSI are very similar to to those derived for the n = CPT-odd 
scenario with the identification in Eq. (JTTJ), (A/^nsi) ~ 7 x 10~ 22 JF GeV. 

More quantitative conclusions on the robustness of the derived ranges for the oscillation 
parameters and on the bounds on the NP strength can be obtained from Figs. 0] and El where 
we plot the marginalized A% 2 as a function of the oscillation parameters, Am 2 and sin 2 8, 
and of the NP strength parameters, for different NP scenarios as labeled in the figures. 

In Table |U we list the 3a allowed ranges for Am 2 and sin 2 9. We read that the derived 
ranges are robust under the presence of these generic forms of NP whose only effect is slightly 
enlarging the allowed range of Am 2 by < 15%, and the lower bound on sin 2 29 by < 7% at 
the 3<r level. We have verified that these conclusions hold for other scenarios characterized 
by different values of n. 

In terms of specific forms of NP the bounds on A5 n imply that ATM+K2K limit the 
possible VLI in the v^—Vr sector via CPT-even effects to 



H < 8.1 x 1(T 25 (1.6 x 1(T 24 ) 



(20) 



and the possible VEP is constrained to 



|0 A 7 | < 4.0 x 10" 25 (8.0 x 10~ 25 ) 



(21) 
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FIG. 4: Dependence of A% 2 on the oscillation parameters Am 2 , sin 2 6 and on the NP strength 
parameter A5 n for different NP scenarios. The full line corresponds to pure —* i/ T Am 2 -OSC. The 
dashed, dotted and dot-dashed lines correspond to different Am 2 -OSC+NP scenarios as labeled in 
the figure. The figure is shown for cr+ = +1. As described in the previous section the results hold 
for <r+ = — 1 with the exchange sin 2 9 — > 1 — sin 2 6 (see discussion around Eq. (|15|0 . The individual 
3d bounds in Tabled can be read from the corresponding panel with the condition A% 2 < 9. 



at 90% (3cr), improving by a factor 8 the previous derived limits in these scenarios [20[. We 
also find that in the v il —v T the VLI via CP-odd effects is constrained to 

\6b\ < 3.2 x lCT 23 (5.0 x 10~ 23 ) GeV (22) 

at 90% (3a), while non-universality of the neutrino couplings to a torsion field 

\Q5k\ < 4.0 x 10" 23 (6.3 x 10~ 23 ) GeV (23) 

at 90% (3a). 

Finally for the case of non-standard neutrino interactions we find the 90% (3<r) bounds 



;-0.021) -0.013 < e < 0.010 (0.017) 
|e| < 0.013 (0.021) 



< 0.029 (0.052) 

< 0.034 (0.060) 



(H NS i real), 
(H NSI complex) 



(24) 



where the upper limits correspond to the case of real NSI and the lower ones to the general 
case of complex e. These limits complement and update the previously derived bounds in 
Refs. 



}f compl 

EQ. 



IV. CONCLUSIONS 

In this work we have studied the robustness of our present determination of the neutrino 
masses and mixing from the analysis of the atmospheric and K2K data under the presence 
of some new physics effects in the v^-Pr sector. In particular, we have performed a global 
analysis to atmospheric and K2K data for scenarios where vector-like or tensor-like new 
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FIG. 5: Dependence of A% 2 on the oscillation parameters Am 2 , sin 2 9 and on the NP strength 
parameters for the case of Am 2 -OSC+NSI. 

physics interactions affect the neutrino evolution together with the standard Am 2 -mixing 
effect. 

We have concluded that the data does not show any evidence of these new physics effects 
even at the sub-dominant level. As a consequence the derived range of oscillation parameters 
is robust under the presence of those unknown effects. The quantification of this statement 
is shown in Figs. H] and 03 and in Table |1J from which we read that inclusion of these new 
physics effects can at most enlarge the allowed range of Am 2 by < 15% and relax the lower 
bound on sin 2 29 by < 7% at the 3a level. 

From the analysis we have also derived upper bounds on the strength of the new physics 
effects in the z^-f T sector. In particular we show in Eqs. (J2U|) and the bound on the 
possible violation of Lorentz Invariance via CPT-even and CPT-odd effects in the neutrino 
evolution respectively. The constraint on the violation of the equivalence principle (VEP) 
due to non universal coupling of the neutrinos to gravitational potential is given in Eq. (|21|). 
while bounds on non-universal couplings of the neutrino to a torsion field are displayed in 
Eq. (j2HJ). The constraints on non-standard neutrino interactions with matter are shown in 
Eq. fl21 • 
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Am 2 -OSC 


Am 2 -OSC+NP 






n = 1 CPT-even 


n = CPT-even 


n = CPT-odd 


NSI 


Am 2 [1(T 3 eV 2 ] 
sin 2 
AS n [GeV n+1 ] 


1.4-3.6 
0.33-0.67 


1.3-3.6 
0.33-0.68 
< 1.6 x 10~ 24 


1.2-3.7 
0.33-0.71 
< 6.3 x 10~ 23 


1.2-3.6 
0.33-0.68 
< 5.0 x 10~ 23 


1.3-3.6 
0.33-0.67 
T < 0.035 



TABLE I: Individual 3a ranges (with 1 d.o.f.) for the oscillation parameters Am 2 and sin 2 9 for the 
different oscillation plus NP scenarios and 3cr bound on the NP strength parameters. The allowed 
range of sin 2 9 corresponds to a£ = +1. For a£ = —1 the corresponding range would be obtained 
by sin 2 9 — > 1 — sin 2 9. 
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APPENDIX A: STATISTICAL TREATMENT OF ATMOSPHERIC DATA 

We summarize here our updated statistical analysis of the atmospheric data. For conve- 
nience we have adopted the pull method used previously by the SK Collaboration (see for 
36| for details on their latest analysis) and by the Bari group 



instance Refs. |35 
There are however some technical differences which we describe next. 

The basic idea of the pull method consists in parametrizing the systematic errors and 
the theoretical uncertainties in terms of a set of variables called pulls, which are then 
treated on the same footing as the other parameters of the model. The \ 2 function can be 
decomposed into the sum of two parts: 

x 2 ($,£) = xLafaQ + xUs(0, (Ai) 

where u denotes the parameters of the model, Xdata * s the usual term describing the deviation 
of the experimental results from their theoretical predictions, and the extra term Xpuiis 
provides proper penalties to account for deviations of the systematics and the theoretical 
inputs from their nominal value. It is convenient to define the pulls in such a way that 
for each source of systematics or theoretical input i the value £j = corresponds to the 
"expected value" reported by the collaboration or predicted by the theory, and & = ±1 
corresponds to a la deviation. 

For the Super-Kamiokande experiment Xp U iis(£) can be properly written as a positive 
quadratic function of £j. The interpretations of the pulls is particularly transparent if the 
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sources of uncertainties are selected to be independent of each other. In this case the pulls 
are uncorrelated and the expression of xpuiis i s ver y simple: 



Xpulls 



(A2) 



In its original formulation, the set of pulls selected by Super-Kamiokande [35] did not 
verify this condition and a correlation matrix between the selected pulls (provided by SK 
collaboration from their MC simulation) had to be included. In our analysis, however, we 
have identified the dominant independent sources of systematic uncertainties in SK analysis, 
and we use them as the basis for our pulls. We have characterized the theoretical and 
systematic uncertainties in terms of 18 independent sources of error: 4 to parametrize the 
theoretical uncertainties associated to the atmospheric fluxes (which we describe in Sec. lA 1|) . 
6 for the theoretical uncertainties in the interaction cross sections (given in Sec. IA 2)1 and 8 
sources of experimental systematic errors (described in Sec. IA 3)1 . To the point to which the 
comparison is possible, this seems close to the approach followed by Super-Kamiokande in 
their latest analysis |36| . 

The form of xlata depends on the expected distribution of the experimental results. Under 
the standard approximation of Gaussian distribution, we have: 



ata 



&3 = £ 



pcx 



a 



stat 



(A3) 



where R^ (R e ^) is the ratio between the expected (observed) number of events and the 
theoretical Monte Carlo for the case of no oscillations. Note that the dependence of Xdata on 
both the parameters uj and the pulls £ is entirely through -R^ h (w, £). In the pull approach, lj 
and £ play a very similar role, and in principle should be treated in the same way. However, 
for the Super-Kamiokande experiment the bounds on £ implied by Xpuiis are m general 
significantly stronger than those implied by Xdata; an d it is therefore a good approximation 
to retain the dependence of Xdata on £ om y to the lowest orders. This is done by expanding 
-R^ h (u;, £) in powers of up to the first order: 



IS 



£<(£)& 



where 



R 



th. 



<V,0), 
dR^&g) 



(A4) 



It is easy to prove 37( that under the approximation (|A4|) Eq. (|A1|) is mathematically 
equivalent to the usual covariance definition of the x 2 which we employed before in Refs. |32l 



391 ] . Thus the small difference in the results are not due to the different statistical treatment, 
but to differences either in the input parameters or in the updated values used for the 
systematic and theoretical uncertainties. 
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Furthermore within the present precision one can safely neglect the dependence of 7r^ on 
the neutrino parameters uj. With this approximation, we can write: 



X \W) = mm 



' 55 

E 

n=l 
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(A5) 



where we have introduced the function x 2 



min{^} X 2 i^i 0- ^ i s c l ear from Eq. (|A5|I 



that in the present approach the systematic and theoretical uncertainties are completely 
characterized by the set of quantities {vr^}, which describe the strength of the "coupling" 
between the pull & and the observable R^. 

In the rest of this section we will discuss in detail how we have parametrized and taken 
into account the various sources of uncertainty and list the derived values for {7i l n }. 



Flux uncertainties 



Flux uncertainties are theoretical uncertainties arising from our limited knowledge of the 



atmospheric neutrino fluxes. Following Refs. 



351 l36| we have parametrized them in terms 



of four pulls: ^j 1 *, £ratio an d ^zenith- The corresponding coefficients ir l n are listed in 

Table HH 



tflux 
>norm 
£j-flux . 



is the pull associated to the total normalization error, which we set to 
=20% 0. 



is a "tilt" factor which parametrizes possible deviations of the energy dependence 
of the atmospheric fluxes from the simple power law. Following Refs 
define: 



?y dep 

3 . 



36 



we 



ME) 



ME) [ I 



ME) 



"4 



(A6) 



and assume an uncertainty on the factor 5, as = 5% Also in analogy with 

Refs. 0,0 we have chosen E = 2 GeV. We then calculate numerically the coeffi- 
cients nll lt as follows: we compute the expected number of events for a given bin N n 
using <&s{E) for the central value of 5 and for 5 ± as and obtain the corresponding 
coupling as the relative change in N n . The results reported in the second column 
of Table ITT1 are obtained neglecting the effect of oscillations. However we have verified 
that when the dependence of the n^ lt on the neutrino oscillation parameters is properly 
taken into account we find very similar results. 

parametrizes the uncertainty on the v^j v e ratio, which is assumed to be a^/ e = 



5% [3fii 0, an d following Ref. we assign a coupling 



-IT, 



C/i/e/2. 



tflux 
t> zenith 



describes the uncertainty on the zenith angle dependence, which we assume 
energy independent. As in Ref. 0] we parametrize the coupling for this pull for the 
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Sample 


Bin 


tflux 
'Snorm 




eflux 
Stilt 






tflux 
>ratio 


tflux 
''zenith 






1 


20% 


(-1.44, 


-1.11) 


x 


5% 

'~J / K) 


(-2 5% 


+2 5%1 


-0.9 x 5 


% 




2 


20% 


(—1.43, 


-1.11) 


x 


5% 


(-2 5% 


+2 5%^ 


—0.7 x 5 


% 




3 


20% 


(—1.42, 


-1.11) 


x 


5% 




+2 5%^ 


—0.5 x 5 


% 




4 


20% 


(—1.42, 


-1.10) 


x 


5% 


(-2 5% 


+2 5%1 


-0.3 x 5 


% 


sub-GeV 


5 


20% 


(—1.42, 


-1.10) 


x 


5% 


(-2.5%, 


+2 5%) 


—0.1 x 5 


% 


(e, u) 


6 


20% 


(-1.42, 


-1.10) 


X 


5% 


(-2.5%, 


+2.5%) 


+0.1 x 5 


% 




7 


20% 


(-1.42, 


-1.10) 


x 


5% 


(-2 5% 


+2 5%) 


+0.3 x 5 


% 




8 


20% 


(-1.43, 


-1.10) 


x 


5% 


C-2 5% 


+2 5%1 


+0.5 x 5 


% 




9 


20% 


(-1.44, 


-1.10) 


x 


5% 


(-2 5% 


+2 5%) 


+0.7 x 5 


% 




10 


20% 


(-1.46, 


-1.10) 


x 


5% 


(-2 5% 


+2 5%) 


+0.9 x 5 


% 




1 


20% 


(+0.35, 


+0.91) 


x 


5% 


(-2 5% 


-1-2 5%1 


-0.9 x 5 


% 




2 


20% 


(+0.38, 


+0.92) 


x 


5% 


(-2 5% 


-1-2 5%"l 


-0.7 x 5 


% 




3 


20% 


(+0.42, 


+0.94) 


x 


5% 


(-2 5% 


+2 5%1 


-0.5 x 5 


% 




4 


20% 


(-1-0.49, 


+0.98) 


x 


5% 


(-2 5% 


+2 5%) 


-0.3 x 5 


% 


multi-GeV 


5 


20% 


(+0.56, 


+1.04) 


x 


5% 


(-2.5%, 


+2 5%) 


-0.1 x 5 


% 


fe. u) 


6 


20% 


(+0.56, 


+1.04) 


x 


5% 


(-2 5% 


+2 5%) 


+0.1 x 5 


% 




7 


20% 


(+0.49, 


+0.98) 


X 


5% 


(-2.5%, 


+2.5%) 


+0.3 x 5 


% 




8 


20% 


(+0.43, 


+0.95) 


X 


5% 


(-2 5% 


+2 5%1 


+0.5 x 5 


% 




9 


20% 


(+0.39, 


+0.93) 


X 


5% 


C-2 5% 


+2 5%1 


+0.7 x 5 


% 




10 


20% 


(+0.35, 


+0.90) 


X 


5% 


(-2.5%, 


+2.5%) 


+0.9 x 5 


% 




1 


20% 




+1.75 x 5% 






-0.9 x 5 


% 


stonninEf 


2 


20% 




+1.72 x 


5% 






-0.7 x 5 


% 


// PVPTlts 


3 


20% 




+1.73 x 


5% 






-0.5 x 5 


% 




4 


20% 




+1.76 x 


5% 






-0.3 x 5 


% 




5 


20% 




+1.84 x 


5% 






-0.1 x 5 


% 




1 


20% 




+4.64 


X 


5% 






-0 95 x 5 


% 




2 


20% 




+4.34 


X 


5% 






—0 85 x 5 


% 




3 


20% 




+4.48 


X 


5% 






—0.75 x 5 


% 




4 


20% 




+4.43 


X 


5% 






—0 65 x 5 


% 


t n m c^m n 


5 


20% 




+4.68 


X 


5% 






-0 55 x 5 


% 


// events 


6 


20% 




+4.62 


X 


5% 






-0.45 x 5 


% 




7 


20% 




+4.61 


X 


5% 






-0.35 x 5 


% 




8 


20% 




+4.96 


X 


5% 






-0.25 x 5 


% 




9 


20% 




+5.01 


X 


5% 






-0.15 x 5 


% 




10 


20% 




+5.22 


X 


5% 






-0.05 x 5 


% 



TABLE II: Coupling factors < of the flux pulls ££ r x m , C r fl atfo and Cnlth with the various 

observables. When the notation (i>i, ^2) is used (second and third column for contained events), 
the first number refer to edike events and the second to //-like events. 
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Sample 


snorm S ra tio 


tl7r tin 
Snorm ?ratio 


tDIS tDIS 
Snorm S ra tio 


sub-GeV e 
sub-GeV \x 
multi-GeV e 
multi-GeV \x 
stopping fj, 
thrugoing ji 


11.3% -0.19% 
11.3% +0.19% 

6.1% -0.20% 

2.1% +0.07% 

2.3% 

0.5% 


3.2% -0.10% 
3.2% +0.11% 
5.0% -0.13% 
5.2% +0.14% 
1.4% 
0.2% 


0.5% -0.01% 
0.5% +0.01% 
3.9% -0.49% 
7.7% +0.98% 
7.5% 
9.6% 



TABLE III: Coupling factors < of the cross-section pulls £? rm, Qtk,, £ no "rm> s^io* Cm and £SI 
with the various observables. The couplings are the same for all the bins in a given data sample. 



bin n of a given sample as ir^ mth = 5% (cos#) n . This means that this uncertainty 
can induce an error in the up/down asymmetry of events which we conservatively take 
to be 5%. In Ref. the assumed up/down uncertainty was smaller (2.5%) and a 
separate zenith-pull was introduced for the horizontal-to-vertical ratio uncertainty of 
2%. We have verified that within the present precision both parametrizations of the 
uncertainties in the zenith angle distribution give very similar results. 



2. Cross-section uncertainties 

Cross section uncertainties are theoretical uncertainties associated to our ignorance on 
the neutrino interaction cross section. In our calculation we follow the standard approach 
and consider separately the contributions to the interaction cross section from the exclusive 
channels of lower multiplicity: quasi-elastic scattering (QE), and single pion production (l7r), 
and include all additional channels as part of the deep inelastic (DIS) cross section (also 
refer to as multi-pion). We neglect for simplicity coherent scattering on oxygen and neutral- 
current interactions, which contribute only marginally to the considered data samples. 

We assume that each of these three contributions to the cross sections are subject to differ- 
ent sources of uncertainties which allow to consider the corresponding pulls as independent. 
For each type of neutrino interactions we introduce two pulls: 

• s^S^m' £norm5 sWm describe the total normalization errors. We conservatively assume 
Cnorm = 15% and cr^orm = 15%. For the normalization error of the DIS cross section 
we estimate cr^orm = 15% for contained events and u^rm = 10% f° r upward-going 
muons from the spread of theoretical predictions arising from the use of different sets 
of nucleon structure functions. The relevant coefficients n l n are listed in Table IIHI 
They are obtained computing the relative change in the number of expected events in 
a given data sample arising from the use of either or ± a^orm f° r eacn °f the three 
contributions to the cross section. 
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Sample 


*sys 
''hadron 


£ S ys 
V/e 


''ring 


?f-vol 


?E-cal 


t s y s 

''FC/PC 


^track 


tsys 
^up-cfT 


sub-GeV e 


-0.25% 


-1.1% 


-0.75% 


-0.3% 


-0.4 


% 







— 


sub-GeV fj, 


+0.25% 


+ 1.1% 


+0.75% 


+0.3% 


+0.4 


% 


— 





— 


multi-GeV e 


-0.50% 


-1.6% 


-2.75% 


-0.5% 


-0.4 


% 








multi-GeV \i 


+1.10% 


+ 1.6% 


+5.40% 


+1.4% 


+2.0 


% 


2.85% 






stopping \i 






+0.30% 


+0.7% 


+0.3 


% 




6.4% 


1.4% 


thrugoing fi 






+0.30% 


+0.7% 


+0.3 


% 




1.4% 


1.0% 



TABLE IV: Coupling factors of the systematics pulls ^ Ton , Q* e , g£g, £f-vol> £&Li» £pc-nrm> 



^sys £ s y s j 

^FC/PC Mrack ana ^up-eff 



tSys f with the various observables. The coefficients are the same for all the bins 
in a given data sample. 



£ratlo' £ratio5 £ratio parametrize the uncertainty of the cri^/a^ ratios. This error is 
relevant only for contained events, and it is much smaller than the total normalization 



uncertainty. The numbers listed in Table IIHI are obtained from Ref. |3_ 



3. Systematic uncertainties 

The systematics uncertainties of the Super-Kamiokande experiment are derived from 
Tables 9.2, 9.3, 9.4 and 9.5 of Ref. 35]. We include in our calculations the following sources 
of systematics: 

Chadron ^ s the P un f° r the uncertainty associated with the simulation of hadronic inter- 
actions; 

£^ is the pull for the errors in the particle identification procedure; 
£ is the pull for the uncertainty coming from the ring-counting procedure; 
£f^oi i s the pull for the uncertainty in the fiducial volume determination; 
C^cai ^ s the pull for the uncertainty in the energy calibration; 

£pc-nrm ^ s the pull for the relative normalization between partially-contained and fully- 
contained events. 

£track ^ s the pull for the uncertainty in the track reconstruction of upgoing muons; 

£up-cfT i s the pull for the uncertainty in the detection efficiency of upgoing muons and 
the stopping-thrugoing separation. 
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